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Abstract. 

We have studied occurrence of quantum phase transition in the one-dimensional 
spin-1/2 Ising model with added Dzyaloshinsky-Moriya (DM) interaction from bi- 
partite and multi-partite entanglement point of view. Using exact numerical solutions, 
we are able to study such systems up to 24 qubits. The minimum of the entanglement 
ratio R = T2/T1 < 1, as a novel estimator of QPT, has been used to detect QPT and 
our calculations have shown that its minimum took place at the critical point. We 
have also shown both the global-entanglement (GE) and multipartite entanglement 
(ME) are maximal at the critical point for the Ising chain with added DM interaction. 
Using matrix product state approach, we have calculated the tangle and concurrence 
of the model and it is able to capture and confirm our numerical experiment result. 
Lack of inversion symmetry in the presence of DM interaction stimulated us to study 
entanglement of three qubits in symmetric and antisymmetric way which brings some 
surprising results. 
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1. INTRODUCTION 

In the last few years it has become apparent that quantum information may lead to 
further insight into other areas of physics such as condensed matter and statistical 
mechanics [Ill2l|3l|4lEll6ll7l|8l|9lll0lllIlll2]. The attention of the quantum information 
community to study in condensed matter has stimulated an exciting cross fertilization 
between the two areas. It has been found that entanglement plays a crucial role in 
the low-temperature physics of many of these systems, particularly in their ground 
statejH m m [71 Ho]. Quantum phase transition (QPT) happens at zero temperature 
and shown non analyticity in the physical properties of the ground state by the change 
of a parameter A of the Hamiltonian H{X). This change is driven only by quantum 
fluctuations ^13j. Since QPT occurs at T = 0, the emerging correlations have a purely 
quantum origin. Therefore, it is reasonable to conjecture that entanglement is a crucial 
ingredient for the occurrence of the QPT |11 El 13 HH [151 US]- Wu et al. [7j have shown 
that a discontinuity in a bipartite entanglement measure (concurrence |15] and negativity 
|16] ) is a necessary and sufficient indicator of a first-order quantum phase transition, 
negativity being characterized by a discontinuity in the first derivative of the ground 
state energy. They have also shown that a discontinuity or a divergence in the first 
derivative of the same measure (assuming it is continuous) is a necessary and sufficient 
indicator of a second-order QPT, that is characterized by a discontinuity or a divergence 
of the second derivative of the ground state energy. 

Dzyaloshinsky has shwon[T7] that, in crystal with no inversion center, the usual 
isotropic exchange JSi.Sj is not the only magnetic interaction and antisymmetric 
exchange Dij.{Si x 5*^) is allowed. Later, Moriya has shown that inclusion of spin 
orbit coupling on magnetic ions in 1st and 2nd order leads to antisymmetric and 
anisotropic exchange respectively. This interaction is, however, rather difficult to 
handle analytically, but it is one of the agents responsible for magnetic frustration. 
Since this interaction may induce spiral spin arrangements in the ground state it is 
closely involved with ferroelectricity in multiferroic spin chains pO|[2Tj. Besides, the DM 
interaction plays an important role in explaining the electron spin resonance experiments 
in some one-dimensional antiferromagnets|22] . Moreover, the DM interaction modifies 
the dynamic properties [23] and quantum entanglement [2^ of spin chains[25]. In the 
present paper, we are interested to study the one- dimensional spin-1/2 Ising model 
with added DM interaction from quantum entanglement point of view using variational 
matrix product state and numerical exact diagonalization methods. The Hamiltonian 
is given by 

N N 

H= jj:s]s^^, +j:d. {s, X , (1) 

where '^j is spin-1/2 operator on the j-th site, and J > (J < 0) denotes 
ant if erro magnetic (ferromagnetic) coupling constant. In very recent works [261 [27] 
respectively studied the ground state phase diagram of the ferromagnetic and 
antiferromagnetic Ising chain in the presence of the uniform DM interaction. It is found 
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that the ground state phase diagram of both systems consists of spiral-ferromagnetic 
and spiral-antiferromagnetic phases respectively and a commensurate-incommensurate 
(C-IC) quantum phase transition occurs at Dc = \ J\. However at the critical value Dc, 
a metamagnetic phase transition occurs into the chiral gapless phase in the ground state 
phase diagram of the ferromagnetic chain. 

This paper is structured as follows. In section 11, we will discuss about bipartite 
and multipartite entanglement measures as QPT indicators of our model and we will 
present our numerical study. In section III, the variational matrix product states will 
be outlined and bipartite entanglement will be obtained. In section IV, we willstudy 
entanglement of three qubits in symmetric and antisymmetric way. Finally we conclude 
and summarize our results in section V. 

2. QUANTUM PHASE TRANSITION 

2.1. Bipartite Entanglement 

The occurrence of collective behavior in many-body quantum systems is associated 
with classical and quantum correlation. The quantum correlation, which known as 
entanglement, cannot be measured in terms of classical physics and represents the 
impossibility of giving a local description of many-body quantum state. The issue of 
finding entanglement measures has recently attracted an increasing interest [Sj El EH |29] . 
Concurrence and tangle are the most widely used measures in QPT related entanglement 
studies. Both of these measures are for bipartite states and because of monogamous 
nature of entanglement they are expected to decrease at the quantum critical point, 
if entanglement is shared by the whole system. Therefore, in order to manifest the 
presence of (QPT) in the model described by Eq.([l]), we focus on the entanglement of 
formation [30] in the quantum spin system and make use of the one-tangle and of the 
concurrence. The one-tangle [HI |32] quantifies the zero temperature entanglement of a 
single spin with the rest of the system and defines as 

ri = 4detp«, p« = i(/ + ^M"^"), (2) 

where p^^^ is the one-site reduced density matrix, = (5*°), and a = x,y, z. In terms 
of the spin expectation values M"', one has: 

ri = l-45:(M"f. (3) 

a 

On other hand, the concurrence [1 5 j quantifies instead the pairwise entanglement 
between two spins and defines as 

Q^ = 2max{0,Ci^,CiS}, (4) 

where 

'-'Im — V ^^^"i yim) + {yim + yim) 
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Figure 1. (Color online.) The concurrence is plotted as a function of DM vector D, 
(a) ferromagnetic and (b) antiferromagnetic cases for different lengths N=12, 16, 20, 
24. 



-\Ji^-9!n.)'-i ' ^ 

ctl = ^{gfZ + 9Z? + {aZ - & 

-^(^.f^)-(^^;^ (5) 

and = (S^S^) is the correlation function between spins on sites / and m and 
Mf = {Sf). The notation (...) represents the ground state expectation value. 
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One-tangle and concurrence are related by Coffman-Kundu-Wootters (CKW) 
conjecture [32], which had been proved by Osborne and Verstraete [33], stating that 

n>r,^Y: CL (6) 

Which expresses the crucial fact that pairwise entanglement does not exhaust the global 
entanglement of the system, as entanglement can also be stored in 3-spin correlations, 
4-spin correlations, and so on. Authors in Ref . ( [28| [29]). have proposed that, due to 
CKW conjecture, the minimum of the entanglement ratio R = T2/T1 < 1, as a. novel 
estimator of QPT, fully based on entanglement quantifiers. 

By doing an experiment, one can find a clear picture of the entanglement 
phenomenon in the ground state magnetic phases of the model. Since a real experiment 
cannot be done at zero temperature, the best way is doing a virtual numerical 
experiment. A very famous and accurate method in the field of the numerical 
experiments is known as the Lanczos method. However, the strong role of a numerical 
experiment to examine quantum phase transitions is not negligible. To explore the 
nature of the entanglement in different magnetic phases, we used Lanczos method to 
diagonalize numerically chains with length up to iV = 24 and coupling constant \J\ = 1. 
The ground state eigenvector, \Gs), was obtained for chains with periodic boundary 
conditions. The numerical Lanczos results on the concurrence for the Ising chain with 
DM interaction, are shown in Fig. [H As is clearly seen from Fig. [Tj^a), in the case 
of ferromagnetic chain and for D < Dc = \ J\ = 1, the concurrence is equal to zero 
which shows that the ground state of the system is in the fully non-entangled polarized 
ferromagnetic phase. At the critical value Dc = \ J\ = 1, the concurrence jumps to a 
non-zero value which confirms the metamagnetic phase transition. By more increasing 
the DM vector, D > Dc, the ground state is in the chiral phase and nearest neighbors are 
entangled. On the other hand, in the case of the antiferromagnetic Ising model, as can 
be clearly seen from Fig. [Tj^b), in the absence of the DM interaction, the ground state is 
non-entangled which is related to the saturated Neel phase. As soon as the DM vector 
applies, nearest neighbors will be entangled and concurrence between them increases 
from zero. Thus in the case of antiferromagnetic Ising chains, the DM interaction induces 
the quantum correlations of the two spins and nearest neighbor spins will be entangled 
as soon as the DM interaction applies. In contrast, in the case of the ferromagnetic 
Ising chains the DM interaction only induces the quantum correlations after the critical 
value Dc and nearest neighbor spins will not be entangled up to the critical value Dc- 

An additional insight into the nature of different phases can be obtained by studying 
the entanglement ratio. Therefore, we have calculated the entanglement ratio by using 
Lanczos method in both ferromagnetic and antiferromagnetic cases. We have plotted 
our numerical results in Fig. [2l For ferromagnetic it can be seen from Fig. |2]^a), 

the entanglement ratio remains zero up to the critical DM interaction Dc which is 
expected from the saturated ferromagnetic phase. As soon as the DM interaction 
increases from the critical Dc, the entanglement ratio starts to increase from zero. In the 
inset of Fig. El^a), the first derivative entanglement ratio is plotted. As it is seen in the 
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Figure 2. (Color online.) The entanglement ratio ti/t2 = R is plotted as a function 
of DM vector D, (a) ferromagnetic and (b) antiferromagnetic cases for different lengths 
N=20, 24. Inset: the first derivative of entanglement ratio. 



ferromagnetic phase, D < D^, the derivative is equal to zero and an abrupt change took 
place exactly at = |^| = 1 which is an indication of the quantum phase transition. 
In the antiferromagnetic case, Fig. EJ^b), as soon as the DM interaction applies the 
entanglement ratio creates and decreases by increasing the DM vector up to = 1-0 
which a change took place and in the D > region the ratio becames monotonous. In 
the inset of Fig. Et^b), the first derivative entanglement ratio is plotted. It is clear that 
in the D > Dc region derivative of ratio is equal to zero and an abrupt change took 
place exactly at Dc = \ J\ = 1.0 which is an indication of the quantum phase transition. 
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2.2. Global Entanglement 



Because of many different kinds of entanglement, quantifying of multipartite 
entanglement states (MES) is more difficult. Global-entanglement {Egi) measure defined 
by Meyer and Wallach [31] which can measure the total nonlocal information per particle 
in a general multipartite system [35] 



gi 



1 

N 



n---iN 



il<i2 



Egi is the average of tangles per particles 



tl<---<lN 
<T> 



(7) 



N 



, without giving detailed knowledge of 



tangle distribution among the individual particles. Therefore, Egi is an average quantity 
and cannot distinguish between entangled states which have equal < r > yet different 
distributions of tangles, like GHZn (Greenberger-Horne-Zeilinger) and EPR®^/"^ state. 
Egi has ability to discriminate between GHZ from W states because of their different 
values of tangle. De Oliveira et al. [36] also introduced a slight extension of global 
entanglement as generalized global entanglement (GGE) which, in contrast to global 
entanglement, the proposed GGE measure can distinguish three paradigmatic entangled 
GHZn (Greenberger-Horne-Zeilinger), EPR®^/'^ and W states. 



N-l 



G(2,n) 



3iV- 



4 
3 



1 



n=l 

3 



N 



N 



j+n 



E (^^^ 

a,/3=0 



13 



(8) 



As such the generalized measure can detect a genuine multipartite entanglement 
and is maximal at the critical point [36]. Here, we have calculated multipartite 
entanglement Eq.(l8]) by using Lanczos numerical method with periodic boundary 
conditions in antiferromagnetic case. We have plotted our numerical results in Fig J3] for 
antiferromagnetic case. As it can be seen from Figj3](a) the multipartite entanglement 
starts to increase by increasing DM up to the critical DM interaction Dc which a change 
took place exactly at Dc = \ J\ = 1 and then after that the multipartite entanglement 
reaches the saturation value. Our calculation shows that multipartite entanglement is 
maximal around the critical point D^. In order to get better insight into the ability 
of multipartite entanglement as a quantum phase transition toolkit we have plotted 
the first derivative entanglement multipartite entanglement in the inset of Fig. [3](a) . It 
shows divergent behavior at the critical point Dc. We have also plotted /im„_j,ooG'(2, n) 
in Fig. |3]^b). Our calculation shows G{2,n) increases as n — )■ oo at the critical point. 
Here, we have finally calculated global entanglement Egi of model Eq.([l]) and compar it 
with multipartite entanglement. It can be seen from Fig. Hlthat both G{2,n) and Egi 
are maximal at the critical point Dc and their behavior is qualitatively the same. 
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Figure 3. (Color online.) (a) Multipartite entanglement estimator G(2, n) is plotted 
as a function of DM vector D for antiferromagnetic Ising chain with length N=20, 
24. Inset: the first derivative of multipartite entanglement '^'^J^'"'' . (b) G{2,n) for the 
chain size N=24, is plotted in the limitc of n — > oo. 



3. VARIATIONAL MATRIX PRODUCT STATE APPROACH 

The matrix product state is defined as|37l |38] 

=Tr{gi,g2,...,gN), (9) 

where gj = aj\ + bj\ Dj, and aj and bj are probabihty amphtudes for two spin 
configurations at site j. In what follow, we intend to determine the ground state energy 
of ferromagnetic Ising spin system with DM interaction. In this respect, by using the 
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Figure 4. (Color online.) Multipartite entanglement G{2, n) and global entanglement 
Egi are plotted as a function of DM vector Z?, for antiferromagnetic Ising chain with 
length N=24. 



above formalism, the variational energy is obtained by 

where (V'lV') = Qi ® Qi-'-Qn ® Qn = ^jGj and Gj = gj (g) gj = |ajp + here 
Hjk = JSjSl + D.{Sj X Sk) and 5" = gj^S^gj. The minimum of the variational energy 
function corresponds to the ground state energy of the system. Using the normalization 
condition = 1 , we can mapped variational parameter to aj = cos{6j)e'^'^\ 

bj = sm{ej)e'^K Therefor one can obtain S'| = |cos(2^j), S+ = {S^Y = | sin(2^j)e^<^> 
where (pj = pj — 0j and by choosing D = Dz, it is found that: 

Evar = ^ Yji-^ cos(26'j) cos(26'j+i) 
j 

+ Dsm{2ej) sm{29j) sm{(pj - (pj+i)}. (11) 

By minimizing above equation, the ground state energy (Eos) in the ferromagnetic case 
J < shall be obtained. It was shown that the ground state energy has the constant 
value, Egs = ~^^^4^ for D < \ J\ and decreasing linearly with DM interaction for D > \ J\ 
as Egs = Now, buy using the minimized variational parameters of E^ar we are 

able to calculate tangle Eq.(|3]) and concurrence Eq.Q. For J > D one can obtain 

XX _ zz _ yx _ _ n 

&i = - 1/4, (12) 
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Figure 5. (Color online.) The ground state energy (normalized by N) of ferromagnetic 
Ising chain with DM interaction using variational matrix product state. Inset shows 
one-tangle and concurrence aas function of the DM interaction using variational matrix 
product state. 



and = so we have ri = 1 and Cj,j+i = 1/2. For J < D again using the above 
conditions, one can obtain 

XX _ yy _ yx _ xy _ 

9t,^i = 1/4, (13) 
and = 1/2 which give ri = and Qj+i = 0. 



4. THREE-QUBIT ENTANGLEMENT 

In this section we focus on the entanglement of formation of three-qubit in two 
inequivalent, symmetric and non-symmetric pairwise entanglement ways. By labeling 
the 3-qubits as 1,2,3 sequentially. The symmetric reduced density matrix pi3 is defined 
as pi3 = tr2{p), where p is the density matrix of 3-qubits. The non-symmetric reduced 
matrix is pi2 = tr^{p). Before present our results, we briefly review the definition of 
concurrence [TSl [32] . Let pij be density matrix of a pair of qubit i and j . The concurrence 
corresponding to density matrix is defined as 

Cij = max{Xi — A2 — A3 — A4, 0}, (14) 

where the quantities Ai > A2 > A3 > A4 are square roots of the eigenvalue of the 
operator 

Qij = PiMy ® (^y)Pij{(^y ® <^y)- (15) 
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The concurrence Ctj = corresponds to an non-entangled state and Cij = 1 corresponds 
to a maximally entanglement sate. A straightforward calculation gives the eigenstates 
and the eigenvalues of 3-qubit of Eq.([T]). The square roots of the operator Qij for 
symmetric and nonsymmetric cases are presented in sequence. In the symmetric case 

^"13 

A. = 0. A, = A3=(^)', A, = 4(-A_)' (16) 

and for the non-symmetric case 

-^-H^f^ --(m^^ ^-(wr <-> 

where a = {J + \/8D'^ + J'^)/2D. In order to determine the existence of entanglement, 
we have considered ferromagnetic and antiferromagnetic cases. Figj6] and FigJT] show 
symmetric, C13, and non-symmetric, C12, concurrences. In the symmetric way. Fig. Ej it 
is clear that in the absence of DM interaction the antiferromagnetic case is entanglement, 
in contrast to ferromagnetic case which is fully unentangled. In antiferromagnetic case, 
the Ci3 is zero up to Dc = \J\ and starts to increase by increasing DM and reaches a 
saturation value. Our calculation shows a competition between Ising exchange J and the 
DM interaction which entanglement starts to decreasing by increasing Ising exchange. 
In the symmetric way, ferromagnetic Ising chain with added DM interaction does not 
show any entanglement and by increasing DM interaction nothing will not happen. 

In the non-symmetric way. Fig. [3, ferromagnetic Ising chain with added DM 
interaction does not show any entanglement up to Z^c = l-'^l- But for > \J\ 
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Figure 7. (Color online.) The ground state entanglement through non-symmetric 
way ,Ci2, for both ferromagnetic an antiferromagnetic cases as function of DM and J 
exchange. 



entanglement starts to increase by increasing DM interaction and shows a competitive 
behavior between DM and J exchange which entanglement shows decreasing by 
increasing \J\. In antiferromagnetic chain, through non-symmetric case, opposite to 
ferromagnetic case, is fully entangled. Entanglement starts to increasing from zero 
as soon as turn on DM and reaches its saturation value around = \J\. It should 
be mentioned that in antiferromagnetic case, increasing J exchange can enhance the 
amount of entanglement. 

The concept of thermal entanglement was introduced and studied within one- 
dimensional isotropic Heisenberg model[39]. Here we study this kind of entanglement 
within three-qubit Ising chain with added DM interaction through symmetric and 
non-symmetric way. The state of the system at thermal equilibrium is p(T) = 
eyopl—H/kT)/^}, where Q = Tr[exp{—H/kT)] is the partition function and k is the 
Boltezmann's constant. As p(T) represents thermal state, the entanglement in the 
state is called thermal entanglement [39] . The square roots of the operator Qij{T) for 
symmetric and nonsymmetric cases are presented in sequence. For the non-symmetric 
case ^13 (T) 



Ai — As 



A3 — (^12 + Z 



12 



A4 



{Y12 - z 



12 



where 



-^12 — — 



12 



-Pes 



1 

2 ^ + 2 



+ 



a 



+ e 
a? + 2 



+ 



+ 



62 + 2 

(6^ + l)e 
62 + 2 
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n2 



1 



a2 + 2 ^ 62 + 2 



and for the non-symmetric case ^12(7") 



Ai — As 



A3 = (ri3 + z 



13) 



A4 = {Y, 



13 



^13j 



(19) 



(20) 



where 



X 



13 



1^13 =^ 



^13 — ~ 



a2 + 2 
(2e- 



1 + 



a2 + 2 



+ 



(2e 



62 + 2 
2e-^^^' 
62 + 2 
2e-^=7' 



+ 2 62 + 2 



(21) 



and n = 4e-^-^[cosh/3J + cosh/3V8£'2 + j2]^ ^ = j _ ^/sD^TJ^/2D, b = J + 
V8-D2 _|_ J2^2D. In FiglH] and Figj9] we give two plots of the thermal concurrence of 
antiferromagnetic and ferromagnetic Ising chains as functions of temperature and DM 
interaction respectively. In antiferromagnetic case, FiglHl model dose not shows any 
entanglement through symmetric way and model lives in fully unentangled phase. In 
contrast to symmetric way, non-symmetric way shows entanglement. As it can be seen 
from FigEl there is a region which surprisingly temperature can enhance entanglement 
and this region will become wide spread for higher DM. For higher temperature DM 
can not overcome temperature and entanglement gradually goes to zero by increasing 
temperature. 

In ferromagnetic case, FiglHl both symmetric and non-symmetric way have 
qualitatively the same behavior. Here, same as antiferromagnetic case, temperature can 
enhance entanglement. Thermal and quantum fluctuations have a tight competition to 
drive system in favorable regime, but in low temperature, both temperature and DM 
interaction has effective influence on the degree of entanglement of the system and both 
of them could be used to increase the entanglement of the spin system. 
Either Figj6] and FigJT] are plotted for D > Dc = I, because we could not find any 
prominent behavior for D < Dc = 1. Actually both symmetric and non-symmetric 
entanglement are zero in the D < = 1 region. 



5. CONCLUSION 



To summarize, we have investigated the effect of a Dzyaloshinskii-Moriya (DM) 
interaction on the ground state phase diagram of the one-dimensional (ID) Ising spin- 
1/2 model using the variational matrix product state and numerical Lanczos methods 
from entanglement point of view. Our results show that there is a critical point in 
either ferromagnetic and antiferromagnetic cases. In the ferromagnetic Ising chain this 
critical point was predicted by variational matrix approach exactly at Dc = \ J\, which 
by using numerical study we have confirmed it. We have also used the minimum of 
the entanglement ratio R = T2/T1 < l,to check the presence of this quantum critical 
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Figure 8. (Color online.) The thermal entanglement through non-symmetric way 
,Ci2, for antiferromagnetic case as function of DM and temperature. 



point. For both ferromagnetic and antiferromagnetic cases our numerical study gave 
the minimum of i? = T2/T1 at Dc = \J\,. We also used generalized multipartite 
entanglement tools to check how entanglement is share in our model and check their 
ability to detect critical point. For antiferromagnetic case we have calculated global 
Egi and generalized multipartite entanglement G{2,n) and either of them show that a 
quantum phase transition took place at Dc = \J\ . 

We have calculated ground state entanglement of symmetric, C13, and non- 
symmetric, C12, concurrences for ferromagnetic and antiferromagnetic cases. In 
the symmetric way, in the absence of DM interaction both antiferromagnetic and 
ferromagnetic cases are fully unentangled. By increasing DM, in antiferromagnetic case, 
the C13 is zero up to critical points Dc = \J\. After the critical point entanglement starts 
to increase until its saturation. For ferromagnetic case, in symmetric way C13 is always 
zero and does not show any entanglement. 

In the non-symmetric way, in the absence of DM interaction, C12 is equal zero 
for either antiferromagnetic and ferromagnetic cases. By turning DM nothing, in the 
ferromagnetic case, will not happen up to critical point Dc and then after that Cu 
starts to increase by increasing DM interaction. In contrast to ferromagnetic case, in 
the antiferromagnetic case Cu starts to increasing immediately after turning DM and 
reaches its saturation value around Dc= \J\. 

We have also studied thermal entanglement of symmetric, Ci3(T), and non- 




Figure 9. The thermal entanglement through non-symmetric way(top) ,Ci2, 
and symmetric way(bottom) ,Ci3,for ferromagnetic case as function of DM and 
temperature. 



J. Vahedi et all. 



16 



symmetric, CuiT) for both ferromagnetic and antiferromagnetic cases. Our calculations 
show that for cither symmetric and non-symmetric cases thermal(uncntanglement 
favorable) and quantum(entanglement favorable) fluctuations have competition to drive 
system to their favorable regime and surprisingly thermal decoherence can enhance 
entanglement in some part of low temperature region. 
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